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ABSTRACT 



Context. Recently, a relationship between the water maser detection rate and far infrared (FIR) flux density has been found as a result 
of a 22 GHz maser survey in a sample comprised of northern galaxies with 100 /jm flux density > 50 Jy and a declination > -30° . 
Aims. The survey has been extended toward galaxies with lower FIR flux densities in order to confirm this correlation and to discover 
additional maser sources for relevant follow-up interferometric studies. 

Methods. A sample of 41 galaxies with 30 Jy < S loo^m < 50 Jy and 6 > -30° was observed with the 100-m telescope at Effelsberg in 
a search for the 22 GHz water vapor line. The average 3<x noise level of the survey is 40 mjy for a 1 km s channel, corresponding to 
a detection threshold for the isotropic maser luminosity of ~ 0.5 Lq at a distance of 25 Mpc. 

Results. Two detections are reported: a megamaser with an isotropic luminosity, Lh 2 o, of ~ 35 Lq in the Seyfert/Hn galaxy NGC 613 
and a kilomaser with L H ,o ~ 1 Lq in the merger system NGC 520. The high luminosity and the presence of a Seyfert nucleus favor 
an association for NGC 613 with an active galactic nucleus. The kilomaser in NGC 520 was also detected with the Very Large Array, 
providing a position with subarcsecond accuracy. The H 2 emission, originating from a £i 0.02 pc sized region with a brightness 
temperature ^ 10 10 K (if the observed variations are intrinsic to the masing cloud(s)), is close to one of the two radio continuum 
sources located in the inner parsecs of NGC 520. The maser is most likely associated with a young supernova remnant (SNR), although 
an association with a low-luminosity AGN (LLAGN) cannot be ruled out. The maser detection rate, with 2 new maser sources out 
of 41 galaxies observed, is consistent with expectations extrapolated from the statistical properties of the S ioo^m > 50 Jy sample. The 
H2O kilomasers are "subluminous", while H9O megamasers tend to be "superluminous" with respect to the FIR luminosity of their 
parent galaxy, when compared with sites of massive star formation in the Milky Way. 

Key words. Masers - Galaxies: active - Galaxies: nuclei - Galaxies: starburst - Galaxies: statistics - Radio lines: galaxies 



1. Introduction 

Water maser emission at 22 GHz (A ~ 1.3 cm) originates from the 
transition between the 616 and 523 rotational levels of ortho-F^O 
and forms in de nse (n(H?) 10 7 cm" 3 ) and warm (7\i n ^;400K) 
molecular gas dHenkel et alj l2005al) . The most plausible in- 
version mechanism is collisional pumping with H2 molecules 
dNeufeld et al.l Il994h . In the Galaxy, water maser sources are 
common in regions of star formation and in the circumstellar en- 
velopes of mass-losing oxygen-rich evolved stars. Since 1977, 
when the first extr agalactic maser was de tected in the nearby 
spiral galaxy M33 ( Churchwell et al.lll977l) . this maser line has 
been observed in ~ 100 external galaxies (Braatz 2007, priv. 
com.). Isotropic 22 GHz line luminosities of extragalactic H7O 
maser s, L H ,o , range from 10~ 2 L q (e. g., M33,[C hurchwell et alj 
119771: IC 10. IHenkel et aDl!986t IC 342. iTarchi et alJl20Q2hh~t o 
23,000 L (SPSS .T0804+3607. iBarvainis & Antonuccil l2005h : 
for a review see[Lo] d2005l) . 

The luminous extragalactic masers, with LhiO > 10 Lq, are 
classified as "megamasers" and are mostly found in galaxies 
with Seyfert 2 or LINER characteristics, exhibiting very high 
X-ray absorbing column densities (TVh > 10 23 ctrT 2 : IZhang et alj 
2006). So far, Very Long Baseline Interferometry (VLBI) studies 
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have been reported for about 10 megamaser sources and all indi- 
cate that the maser emission is located within a few parsecs from 
the associated active galactic nucleus (AG N), being related to (1) 
nuclear accretion dis ks (e.g., NGC 4258,lGreenhill et al.lll995 , 
i Mivoshi et al1[T995t IC 2560. Ilshihara et alj l200lt Mrkl419, 
iHenkel et alj|2002l) . (2) interactions between the jet(s) emerg- 
ing from t he nucleus and amb ient molecula r clouds (e.g., 
NGC 1052, IClaussenet all II 998: NGC 1068, iGallimore et ail 
l2Q0lh Mr k348. iPecket alJ 120031) . and (3) nuclear outflows 
(Circinus, iGreenhill et alj|2003al) . 

The less luminous extragalactic masers, the so-called 
"kilomasers", with Lh,o<10Lq, are usually found in star- 
forming galaxies and the majority of them are associated 
with star-for ming regions clearly offset from the nucleus 
(e. g. M 33, IChurchwell etafj|1977t M 82 iBaudrv & Brouilletl 
1996; IC1 IHenkel et alj fl986r. IC 342. ITarchi et all l2002bt 
NGC 2146. [Tarchi et al.l l 2002al) . H owever, in 3 galaxies M 5 1 
( Hag iwara et all 1200 ll) NGC 405 1 dHagiwaraet alj|2003l) . and 
NGC 253 dHenkel et alj|2004l) . kilomaser emission is also aris- 
ing from the nuclear region, with projected distances from the 
putative nucleus of < 15, 5, and 6 pc, respectively. Among these, 
only the kilomaser in M 5 1 has been definitely assoc iated with 
an AGN, most likely with the radio jet dHagiwaral2007l) . Instead, 
the kilomaser in NGC 253 turned out to be associated with a su- 
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perno va remnant in the proximity of the nucleus dHofner et al.l 
2006). The nature of the water maser emission in NGC4051 is 
still uncertain an d needs to be investigated at higher resolution 
(Hagiwara 2007). The possibility that some kilomasers could be 
related to nuclear activity (like M51 and, possibly NGC4051), 
similarly to the more powerful megamasers, is exciting and 
would provide a new tool for studying their host galaxies' nu- 
clear regions. Hence, given the small number (15) of kilomasers 
found so far, this motivates extended searches to find more of 
these sources. 



The sample: Searching for 22 GHz water vapor masers toward 
a sample comprised of all galaxies with declination S > -30° 
and 100;um IRAS point so urce flux densities Si oo um >50Jy, 
led to four new detections (Tare hi et al.l I2002alfb1 : iHenkel et al.l 
l2005bl hereafter HPT), yielding a total detection rate of 22%. 
A relationship between FIR flux density and likelihood of de- 
tectable water maser emission was also established. Motivated 
by these results, we have decided to extend the observational 
campaign toward lower infrared fluxes, to confirm the corre- 
lation found and to discover additional maser sources suitable 
for interferometric follow-up observations. The sample includes 
all galaxies with 6 > -30° and IRAS point source flux den- 
sities 30 Jy< Sioo//m < 50 Jy with known recessional velocities. 
The list of sources, shown in Table[TJ was compiled usin g the 
IRAS Point Source Catalog dFullmer & Lonsdald Il989l). For 
three of the 41 target ga laxies listed, NGC 613 (IHenkel et al.l 
ll984t lBraatzet all (l 9961). NGC 5 135 dBraatz et al.l \ 2Q04. and 
NGC 6951 (IClaussen & Lolll986t lBraatzetal.lll996ll2004l) up- 
per limits had already been determined. Because of the time 
variability of the 22 GHz H2O maser emission, the sources were 
re-observed. In the following, we describe the results of this new 
survey. 



2. Observations and data reduction 

2.1. Effelsberg observations 

Measurements in the 616 - 523 transition of ortho-H20 (rest fre- 
quency 22.23508 GHz) were carried out with the the MPIfR 
100-m telescope at Effelsberg, Germany, between April 2004 
and November 2005. The full width to half power (FWHP) 
beamwidth was w 40" and the pointing accuracy was, on av- 
erage, 4" and always better than 20". Throughout the survey, 
we used a K-band high electron mobility transistor (HEMT) re- 
ceiver with two orthogonally linearly polarized channels and a 
rotating feed that was switched between two fixed position sep- 
arated by ~ 2' at a frequency of 1 Hz. The spectrometer was an 
8192 channel correlator (AK 90) in which the signal was split 
into 8 bands (IFs) of width 40 MHz and 512 channels each, 
yielding a channel spacing of ~ 78 kHz (~ 1.04km s _1 ). Four in- 
termediate frequency bands ("IFs") were centered at a frequency 
corresponding to the local standard of rest (LSR, optical conven- 
tion) recessional velocity of each galaxy, while the others had 
frequency offsets of +20 MHz. On January 21, 2005 we also 
took a high spectral resolution spectrum of NGC 520 using 2 
bands of 20 MHz and 4096 channels each, yielding a spacing of 
~ 5 kHz (0.07 km s _1 ). 

To establish a flux density scale, the continuum sources 
NGC 70 27, 3C 286, and W3 (OH) were observed (for fl ux densi- 
ties, see lOtt et alll 19941 and iMauersberger et al.ll 1988b . All data 



were reduced using standard procedures with the CLASS pro- 
gram of the GILDAS software package. For each source, all 
scans at the same frequency were added to reduce the noise in 
the final spectrum. All spectra were calibrated to convert the in- 
tensity scale from noise diode uni ts to flux density (in Jy ), using 
an updated version of Eq. 1 of Gallimore et alJ d200l[) . which 
also takes into account gain variations as a function of elevation: 



[Jy] 



N 



2.5 



[counts] [Jy counts '] 

{ [1 -0.01(30° -0°)Y 
x I 1.0, 

I [1 -0.01(0° -60°)]" 



x 



if 6° < 30°, 

if 30° < 6° < 60° 

if 6° > 60°, 



where: 



- N is the signal measured by the Effelsberg telescope 

- the third term corrects for the elevation dependent gain of the 
telescope (6°: elevation in degrees). 

The flux calibration uncertainty is estimated to be of the order of 
+ 10%. More observational details are given in TableQ] 



2.2. VLA observations of NGC 520, archival data, and 
imaging 

NGC 520 was observed at 22 GHz with the NRAC0 Very Large 
Array (VLA) in the BnA and A configuration on February 14, 
2005 and April 17, 2006, respectively. In February 2005 obser- 
vations were made in Spectral Line Mode, using a single band- 
width of 3.1 MHz, which corresponds to a velocity coverage of 
~45km s _1 , centered at VLSR=2276km s _1 . The use of 256 
spectral channels yielded a spacing of 12.2 kHz, corresponding 
to 0.16 km s" 1 . In April 2006, NGC 520 was observed both in 
Spectral and Continuum Mode. For the spectral line observa- 
tions, we employed the set up used in February 2005, while 
continuum observations were made using two bands of 50 MHz 
each, centered at frequencies of 22.4851 and 22.4351 GHz. 
"Referenced Pointing" observations were made in regular in- 
tervals to warrant good antenna pointing. The "Fast Switching" 
technique was used with a typical duty cycle of 4 minutes, 
~170sec on the target source (NGC 520) and ~70sec on the 
phase calibrator 0121+043 (0.67Jy). The flux density scale 
was established by observing the non-variable source 0137+331 
(3C48; 1.13 Jy). 

Archival VLA A-array continuum data at 8.5 GHz (project 
AA296, P.I. Hofner) and 14.9 GHz (project AC156, P.I. Carrol) 
were also reduced by us. These observations were taken in 
December 2004 and May 1986, respectively. Absolute calibra- 
tion was established at 8.5 GHz using 3C48 (3.15 Jy) and at 
14.9 GHz using 3C286 (3.45 Jy). The compact sources 0125- 
000 (1.15 Jy) and 0106+013 (3.11 Jy) were measured to cali- 
brate the telescope phases. Furthermore, we have produced two 
continuum MERLIN^ maps at 1.4 and 5.0GHz. NGC 520 was 
observed at 1.4 GHz and 5.0 GHz with the full MERLIN array 
in December 2000 and November 2004, respectively. For the 1 .4 
GHz session, the 76-m Lovell telescope took also part in the ob- 
servations. 



1 The Effelsberg 100-m telescope is a facility of the MPIfR (Max- 
Planck-Institut fur Radioastronomie) in Bonn. 



- The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

3 MERLIN is operated by the University of Manchester as a National 
facility of the Science and Technology Facilities Council (STFC). 
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Table 1. List of galaxies searched for H2O maser emission. 
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Source 




RA^ 




Dec" 


"sys 


c b 


V-range^ 


rms 


La 


d 

:Q_ 


Epoch 




(J2000) 




(J2000) 


(kms^ 1 ) 


(Jy) 


(km s _1 ) 


(mJy) 


L ( 


3 




NGC 157 


00 


34 


46 


.7 


-08 23 47 


1668 


37 


.73 


1120, 


2220 


10 


<0 


.39 


21-Jan-05 


NGC 278 


00 


52 


04 


.3 


47 33 02 


640 


45 


.16 


100, 


1180 


8 


<0 


.04 


21-Jan-05 


NGC 52<K 


01 


24 


35 


.1 


03 47 33 


2281 


48 


.40 


see Table [3] 










NGC 613~ 


01 


34 


18 


.2 


-29 25 07 


1475 


49 


.14 


see Table [3] 










NGC 908 


02 


23 


04 


.6 


-21 14 02 


1509 


44 


.43 


960, 


2060 


13 


<0 


.40 


21-Jan-05 


NGC 1385 


03 


37 


28 


.3 


-24 30 05 


1493 


35 


.62 


950, 


2040 


22 


<0 


.65 


19-Jan-05 


UGC 2866 


03 


50 


14 


.9 


70 05 41 


1232 


49 


.99 


690, 


1780 


14 


<0 


.29 


06-Apr-04 


NGC 1482 


03 


54 


38 


.9 


-20 30 09 


1916 


46 


.52 


1370, 


2470 


19 


<0 


,96 


10-Apr-05 


NGC 1614 


04 


33 


59 


.8 


-08 34 44 


4778 


32 


.38 


4210, 


5340 


14 


<3 


.44 


12-Apr-04 


NGC 2339 


07 


08 


20 


.5 


18 46 49 


2206 


30 


.24 


1660, 


2760 


9 


<0 


.59 


19-Jan-05 


NGC 2566 


08 


18 


45 


.7 


-25 30 03 


1637 


39 


.43 


1130, 


2220 


30 


<1 


.12 


12-Apr-04 


NGC 3175 


10 


14 


42 


.1 


-28 52 19 


1101 


31 


.44 


560, 


1650 


27 


<0 


.44 


21-Jan-05 


NGC 3310 


10 


38 


45 


.9 


53 30 12 


993 


41 


.76 


450, 


1540 


13 


<0 


.17 


06-Apr-04 


MESSIER 95 


10 


43 


57 


.7 


11 42 13 


778 


35 


.30 


230, 


1320 


9 


<0 


.07 


21-Jan-05 


NGC 3504 


1 1 


03 


11 


.2 


27 58 21 


1534 


32 


.89 


990, 


2080 


8 


<0 


.24 


21-Jan-05 


NGC 3593 


1 1 


14 


37 


.0 


12 49 05 


628 


35 


.79 


90, 


1170 


8 


<0 


.04 


21-Jan-05 


NGC 3675 


1 1 


26 


08 


.6 


43 35 09 


770 


35 


.20 


230, 


1310 


6 


<0 


.05 


21-Jan-05 


NGC 3810 


1 1 


40 


58 


.8 


11 28 17 


993 


31 


.40 


450, 


1540 


9 


<0 


.12 


31-Mar-05 


NGC 3893 


1 1 


48 


38 


.2 


48 42 39 


967 


35 


.07 


420, 


1510 


9 


<0 


.11 


31-Mar-05 


NGC 4030 


12 


00 


23 


.6 


-01 06 00 


1465 


46 


.34 


920, 


2010 


10 


<0 


.28 


21-Jan-05 


NGC 4041 


12 


02 


12 


.2 


62 08 14 


1234 


31 


.49 


690, 


1780 


14 


<0 


.29 


06-Apr-04 


NGC 4157 


12 


1 1 


04 


.4 


50 29 05 


774 


43 


.24 


230, 


1320 


8 


<0 


.07 


10-Apr-05 


NGC 4355 


12 


26 


54 


.6 


-00 52 39 


2179 


33 


.78 


1630, 


2730 


12 


<0 


.79 


19-Jan-05 


NGC 4536 


12 


34 


27 


.1 


02 11 16 


1808 


44 


.98 


1260, 


2360 


8 


<0 


.36 


31-Mar-05 


NGC 4605 


12 


39 


59 


.4 


61 36 33 


143 


30 


.53 


-400, 680 


13 


<0.i 


303 


06-Apr-04 


NGC 4654 


12 


43 


56 


.6 


13 07 35 


1037 


35 


.26 


490, 


1580 


11 


<0 


.16 


10-Apr-05 


UGC 8058 


12 


56 


14 


.2 


56 52 25 


12642 


30 


.89 


12030, 13250 


11 


<21 


.57 


06-Apr-04 


NGC 5033 


13 


13 


27 


.5 


36 35 38 


875 


44 


.53 


330, 


1420 


11 


<0 


.11 


19-Jan-05 


NGC 5078 


13 


19 


50 


.0 


-27 24 36 


2168 


32 


.76 


1620, 


2720 


30 


<1 


,93 


12-Apr-04 


NGC 5135 


13 


25 


44 


.0 


-29 50 01 


4112 


30 


.83 


3550, 


4670 


30 


<6 


.22 


ll-Apr-04 


NGC 5248 


13 


37 


32 


.1 


08 53 06 


1153 


43 


.97 


610, 


1700 


10 


<0 


.19 


21-Jan-05 


MESSIER 101 


14 


03 


12 


.6 


54 20 57 


241 


30 


.72 


-300, 780 


17 


<0 


.01 


06-Apr-04 


NGC 5676 


14 


32 


46 


.8 


49 27 28 


2114 


30 


.61 


1560, 


2670 


9 


<0 


.57 


21-Jan-05 


NGC 5713 


14 


40 


11 


.5 


-00 17 21 


1899 


37 


.22 


1350, 


2450 


19 


<0 


.92 


21-Jan-05 


NGC 5775 


14 


53 


57 


.6 


03 32 40 


1681 


45 


.21 


1130, 


2230 


7 


<0 


.26 


31-Mar-05 


NGC 5907 


15 


15 


53 


.7 


56 19 44 


667 


36 


.01 


120, 


1210 


11 


<0 


.07 


07-Apr-04 


NGC 6643 


18 


19 


46 


.4 


74 34 06 


1484 


31 


.78 


940, 


2050 


15 


<0 


.44 


10-Apr-05 


NGC 6951 


20 


37 


14 


.1 


66 06 20 


1424 


37 


.14 


880, 


1970 


12 


<0 


.33 


10-Apr-05 


NGC 7469 


23 


03 


15 


.6 


08 52 26 


4892 


35 


.22 


4330, 


5460 


8 


<2 


.34 


12-Apr-04 


NGC 7541 


23 


14 


43 


.9 


04 32 04 


2678 


39 


.93 


2120, 


3230 


11 


<1 


.05 


12-Apr-04 


NGC 7771 


23 


51 


24 


.9 


20 06 43 


4277 


38 


.78 


3710, 


4840 


11 


<2. 


.47 


12-Apr-04 



" Coordinates and velocities are taken from t he NASA/IPAC Extraealacti c Database (NED). The systemic velocities are given by V sys = cz. The 
exception is NGC 613 for which the velocity of de Vaucoul eurs et al.l il991) has been taken (see also Sect. l3.lt . 
* Flux densities at 100 yum are taken from the IRAS Point Source Catalog (Fullmer & Lonsdale 1989). 
c Inspected velocity ranges. 

d Upper limits for maser luminosities are derived from Lh,o/[L0] = 0.023 x 5/[Jy] x Av/[km s~'] x £> 2 /[Mpc 2 ], where S is 3 times the rms in 
column 7 and Av is the channel width, 78 kHz or ~ 1 .04 km s . 
e Source name in bold-face: newly detected maser. 

Table 2. Interferometric continuum maps of NGC 520. 



Band 


V 


Observation 


Telescope 


Weighting 


Restoring Beam 


Position 


n a 


rms 




(GHz) 


date 






HPBW (") 


Angle O 


(") 


(mJy/fi,,) 


L 


1.4 


10-Dec-2000 


MERLIN 


Natural 


0.350 x 0.350 


0.00 


3 


0.1 


C 


5.0 


19-Nov-2004 


MERLIN 


Tapered 1.3 Ml 


0.187x0.123 


-51.12 


1 


0.07 


X 


8.5 


27-Dec-2004 


VLA-A 


Uniform 


0.223 x 0.215 


18.89 


7 


0.03 


u 


14.9 


10-May-1986 


VLA-A 


Natural 


0.205 x0.151 


19.81 


4 


0.1 


K 


22.1 


14-Feb-2005 


VLA-BnA 


Natural 


0.369x0.136 


63.21 


7 


0.5 


K 


22.5 


17-Apr-2006 


VLA-A 


Natural 


0.112x0.084 


11.29 


1 


0.06 



" Largest angular scales that can be imaged. 
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All datasets were calibrated in the standard way, except that, 
since 3C48 is resolved by the VLA both at 8 and 22 GHz, 
phase and amplitude correction factors were determined by us- 
ing CLEAN component models for this source. Each individ- 
ual dataset was Fourier-transformed using natural weighting 
(with the exception of the 8.5 GHz map where we used uniform 
weighting) and then deconvolved using the CLEAN algorithm 
(Hogbom 19741) . Table [2]illustrates the details of the continuum 
maps. 



3. Results 

Toward two of the 41 observed targets, the Seyfert galaxy 
NGC613 and the merging system NGC520, H2O maser emis- 
sion was detected. Figs.[T]and|2]show the line profiles, while line 
parameters are given in Table[3] Below, we summarize the main 
properties of the two galaxies and their masers. 



3.1. NGC613 



NGC 613 is a barred spiral galaxy dde Vaucouleurs et alj 
199 ll) . located at a distance of ~20Mpc (1 " = lOOpc, 
for i/o=75kmsMpc _1 ). Optical spectroscopy revealed that 
this galaxy has a composite spectrum, with an Hn re- 
ion and a Seyfert-like co ntribution simultaneously present 
Veron-Cettv & Veronlll986l) . High resolution radio continuum 
maps of the inner kpcs of NGC613 show a linear structure, 
consisting of three distinct blobs, oriented perpendicular to a 
ring-like feature of radius ~350pc (Hum mel & Jorsater|[l 992). 
A large amount of CO has been detected dElfhag et al . 1996) but 
high angular resolution images of molecular emission have, so 
far, not yet been reported. 

The H2O maser profile in NGC613, shows two distinct 
components (see Fig. [1). The most luminous one is centered at 
~ 1500km s _1 , which is redshifted w.r.t. the systemic velocity 
of the galaxy, and has a full width to half maximum linewidth 
of ~40km s _1 . A second blueshifted component is seen at 
~ 1450km s _1 and has a linewidth of ~20km s _1 . The total 
isotropic maser luminosity is ~ 35 Lq. Independently from our 
observation, the maser was detected by iKondratko et all {2006) 
in July 2003. They report a position for the maser coincident 
with the optical position of the nucleus with an uncertainty of 
1 "3, supporting the association of the maser emission with AGN 
activ ity. The comparison of our spectrum, from 2005, with that 
of IKondratko et al.1 (120061) (their Fig. 1), from 2003, shows that 
the line profiles are similar. However, in the former, the peak 
flux density of the component at Vlsr ~ 1450km s appears 
to be higher. The large linewidth and the presence of not more 
than two main components, separated by only a few 10 km s" 1 , 
may hint at a jet-maser, as proposed by Kondrat ko et al.l d2006). 
Nevertheless, the H2O emission is centered on the Hi systemic 
velo city (see Fig. JT) an d on the CO centroid velocity (1465 km 
s ; lElfhag et al.l 1 1996b . which is not common for jet-masers. 
Characteristically, jet-masers are entirely found at one side of 
the systemic velocity, either blue- or redshifted 

dEB). 



3.2. NGC 520 

NGC520 is a nearby (D~30Mpc, 1" = 145 pc, for 
//o=75kmsMpc _I ) merging system, considered to be the re- 
sult of an encounter between a gas rich and a gas poor 
disk galaxy dHibbard & van Gorkom|[l996l) . The remnant nu- 
clei are clearly visible in near infrared K-band (2 pm) images 
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Fig. 1. H2O megamaser profile observed toward NGC 613 with 
Effelsberg on January 19, 2005. Channel spacing and pointing 
position are 4.3 km s _I and 02000 = 01 h 34 m 18 s . 2, £2000 = - 
29° 25' 07", respectively. The velocity scale is with respect to 
the local standard of rest (LSR) and uses the optical convention. 
The vertical line indicates the Hi systemic velocity V sys = 1481 + 
5 km s- 1 (NED). 



dStanford & Balcellsl 119901: iKotilainen etaTJ 1200 lb . At optical 
wavelengths, the pri mary nucleus, which is fo ur times more mas- 
sive than the other (Stanford & Balcells 1990), is hidden behind 
a pronounced dust lane, while the s econdary less massive one 
is seen 40" (6kpc ) to the north-west dStanford & Balcells|[l990t 
lLaine et alj|2003h . 

The radio continuum emission arising from NGC 520 is con- 
fined to a linear east- west structure of ~ 6"(~ 870 pc) in extent, 
with the peak coinciding wi t h the 2/j.m position of the main 



witn tne peak coinciding wit n tne z/i m positio n 01 tne marr 
nucle us dCondon etalj|1982t ICarral et al.ll 19901: iBeswick et al 



2003). The 4.9 GHz unresolved peak source of ICondon et al 
d 19821) is resolved into 3 components in the 14.9 GHz map of 
ICarral et ail d!990h . The 14.9 GHz map re veals a total of 5 com- 
ponents, aligned along the linear structure. IBeswick et alj (2003) 
imaged NGC 520 at 1.4 and 1.6 GHz with MERLIN. In their 
highest resolution map, the linear continuum structure resolves 
into 10-15 compact components. Five o f these co incide with the 
source s detected by ICarral et"aT] d 1990b . which IBeswick et al.l 
(I2003h associated with bright supernova remnants (SNRs) and 
radio supernovae (RSNe). 

The inner 1 kpc surrounding the primary nucleus contains a 
large concentration of neutral atomic and molecular gas. Hi and 
OH hav e been detected in ab sorption, along with an OH maser 
feature dBeswick et al.l l2003). The velocity distribution of both 
Hi and OH trace a kpc-scale ring or disc at a position angle (P. A.) 
of 90° rotating about the main nucleus with a velocity gradient 
of 0.5 km s' 1 pc" 1 . This resu lt is consistent with the CO velocity 
field dYun & Hibbardll2001l) . 

Chandra X-ray observations (Read 2005) show no signifi- 
cant diffuse emission from the secondary nucleus, while diffuse 
emission, consistent with a starburst driven wind, has been dis- 
covered outflowing fro m the primary nucleus , perpendicular to 
the molecular CO disk dYun & Hibbardll2001l) . Chandra obser- 
vations also revealed 15 point so urces within the optical bound- 
aries of the galaxy (Read 2005). The third brightest of these, 
source 12, is an ultraluminous X-ray source (ULX) with Lx(0.2- 
10keV)~ 9xl0 39 ergs~ 1 . S ince it coincides, wi thin the uncer- 
tainty, with the radio peak (Condon et al. 1982 ), the center of 
the molecular disk traced by CO dYun & Hib bard 2001), and the 
center of the diffuse wind structure, it has been associated with 
the main nucleus of NGC 520. 
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Source 


Adopted 
Distance^ 
(Mpc)~ 


Epoch 


Tel. 


J SdV^ 
(mJy km s~') 


V a 

v LSR.opt 

(km s" 1 ) 


"'1/2 
(kms- 1 ) 


1 j h 

K-range_ 

(kms^ 1 ) 


Channel 
Width 

(km s" 1 ) 


rms 
(mJy) 


(L Q ) 


NGC 613 


20 


Jan 19 2005 


Eff 


1500±300 


1453+2 


23±4 


930, 2020 


4.3 


18 


14 










2300±300 


1502+3 


37±6 




4.3 


18 


21 


NGC 520 


30 


Jan 21 2005 


Erf 


120±10 


2275.28+0.05 


1.1+0.6 


1730, 2830 


1.1 


8 


2.1 










70±60 


2275.42±0.08 


0.6±0.2 


2010, 2560 


0.3 


28 


1.4 






Feb 14 2005 


VLA 


25 ±2 


2275.15±0.03 


0.61+0.07 


2250, 2300 


0.2 


4 


0.5 






Mar 31 2005 


Eff 








2010, 2560 


0.3 


29 


<0.5 






Apr 10 2005 


Eff 








1730, 2830 


1.1 


11 


<0.8 






Nov 14 2005 


Eff 








2010, 2560 


1.1 


8 


<0.5 






Apr 17 2006 


VLA 








2250, 2300 


0.2 


5 


<0.1 



Integrated flux densities, center velocities, and the full width to half maximum (FWHM) linewidths of the observed features are obtained from 
Gaussian fits. 
* Inspected velocity ranges. 

c Isotropic luminosities are derived from L H W[L©] = 0.023 x f SdV/[Jy km S" 1 ] x D 2 /[Mpc 2 ]. 
d Distances have been estimated using H =75 km s -1 Mpc -1 . 
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Fig. 2. H2O kilomaser spectra observed toward NGC 520 be- 
tween January 21 and March 31, 2005. Channel spacings are 
0.6kms _I (lower panel), 0.2 km s _1 (middle panel), and 1.1 km 
s _1 (upper panel). The pointing position is (22000 = 01 h 24 m 35 s . 1, 
£2000 = +03° 47' 33". The solid vertical line indicates the sys- 
temic velocity V sys =2281km s _I (NED). The dash-dotted and 
dashed lines represent the peak velocities of the maser in the 
Effelsberg and VLA spectra, respectively. 



3.2.1 . The detection and variability of the H2O maser 

On January 21, 2005, we detected a single, narrow H2O maser 
line in NGC 520 (see Fig |2). The line parameters are reported 



in Table [3] The inferred isotropic maser luminosity is ~ 1 L0. 
On February 14, the VLA observations confirmed the presence 
of a single, narrow emission feature. As is evident from Fig. [2] 
and Table [3j the peak flux density of the line has decreased 
w.r.t. the earlier Effelsberg spectrum. In addition, a small shift 
of ~ 0.3 km s toward the blue is apparent. A 2-dimensional 
Gaussian fit (made using AIPS task JMFIT) to an image of 
the channel maps showing maser emission indicates a spa- 
tially unresolved spot at a position a 2 ooo=01 h 24 m 34 s .914±0 s .001 
and £ 2 ooo=03°47'29'.'79+0''01. The uncertainties of the position 
are the formal errors determined by JMFIT. These are consis- 
tent with the theoretical error given, for each coordinate, by 
~0.59b/S NR, where 9b is the beam (HPBW) size along either 
RA or Dec and SNR is the signal-to-noise ratio in the map. 
Having used the 'Fast Switching' technique (see Sect. |2.2t , the 
absolute positional accuracy is likely consistent with this value. 
Nevertheless, for the absolute uncertainty we conservatively take 
O'.'l, the canonical error for VLA K-band A-array observations. 
The position and the velocity of the H2O mas er are incompatible 
with those of the weak OH masers detected bv lBaan et alJ (Tl985) 
and lBeswick et alj d2003l) . 

In the Effelsberg observation made on March 3 1 , no emission 
above the 3c noise level of ~90mJy was found in the 0.3 km 
s wide channels (Fig.|2]shows a smoothed spectrum). The fad- 
ing of the maser line was confirmed by our most recent spec- 
trum, taken with the VLA on April 17, 2006, with a 3cr detection 
threshold of 15 mJy for a 0.2 km s wide channel. The decrease 
in flux density of ~70% within 24 days (from January 21 to 
February 14), implies a size scale $6.2xl0 16 cm (S0.02pc) or 
^ 0.14 mas at a distance of 30 Mpc (if variations are intrinsic to 
the masing cloud(s)). The corresponding brightness temperature 
is>10 10 K. 

We produced continuum maps from the VLA line data cubes 
of February 2005 and April 2006 using the line free channels 
and find an unresolved spot with a peak flux density of 3- 
4mJybeairT 1 at the position of the maser. We compared these 
maps with the broad band continuum map at 22. 1 GHz (see next 
Sect.). Position and flux density of this source are not compat- 
ible with those of the continuum source detected in the broad 
band map. In order to investigate this ambiguity, we smoothed 
the line data cubes to a resolution of 0.8 km s to increase the 
signal-to-noise ratio in individual channels. In these maps we 
see "tentative" weak line emission at a 3-4cr level at exactly the 
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Fig. 3. VLA A-array 22.5 GHz broad band continuum map of 
NGC520, from April 17, 2006. Contour levels are -3, 3, 4, 5, 
6 x 0.06 mJy beam 1 . The cross marks the position of the H2O 
maser. 

same position as the main maser line but at different velocities. 
We made additional continuum maps excluding these velocity 
channels and the continuum source is still present at the same 
position, though with lower flux density. Hence, we think that 
the spot seen in the line-free maps is likely caused by the blend- 
ing of weak maser lines below the detection threshold of the 
individual channel maps, possibly spread over the entire 45 km 
s _1 bandpass. 

3.2.2. The 22 GHz continuum source 

In the 22.5 GHz continuum image made from the broad 
band continuum dataset, we find a single, almost unre- 
solved source (Fig. [3j hereafter labeled "K"), whose posi- 
tion, determined fitting an elliptical Gaussian with JMFIT, is 
ar 2 ooo=01 h 24 m 34^921+0 s :001 and 5 2 ooo=03°47'29'.'73+0'.'01. If 
we assume, as before, that the absolute positional accuracy is 
O'.'l, then the uncertainty on the relative positions is 0'.'14. 
Within this conservatively estimated error the maser and the con- 
tinuum source coincide. However, as we stated in Section [3.2.1l 
the absolute position error is probably smaller then this value 
and might even be of order 10 mas. In this case the peaks of 
the maser and that of source K do not coincide any longer (they 
are displaced by 0'.'12±0"014, or 18+2 pc). The peak and in- 
tegrated flux densities of source K are 0.41+0.07 mJybeanT 1 
and 0.5+0.1 mJy, respectively. The undeconvolved extent is 
110+20xl00+20mas with position angle (P.A.) ~ +30°, match- 
ing the size and P.A. of the beam. The brightness temperature of 
the source is <: 40 K. 

3.2.3. Nuclear continuum sources 

Figure [4] shows the MERLIN and VLA images of the inner 
6"x275 region surrounding the main nucleus of NGC520 at 
the frequencies of 1.4, 5.0, 8.5, and 14.9 GHz (details of these 
maps are given in Tab le [2}. Consistent with p revious studies 
(e.g. ICarraletal.il 19901 anc UBeswick et all 120031) . a linear east- 
west structure is visible at all frequencies. This structure is re- 
solved into several compact components, the number of which 
depends on the resolution and sensitivity of the maps. In our 



maps it is possible to identify a total of 7 compact sources, la- 
beled with the letters A to G. The four sources A, B, C, and F, are 
resolved into sub-components in the MERLIN C-band and VLA 
U-band maps, because of their higher resolution. Coordinates, 
peak and integrated flux densities are displayed in Table [4] for 
the MERLIN maps and in Table [5] for the VLA data. From a 
comparison with previously published radio continuum maps, 
we can associate our five s ources A - E with those lab e led C I 
- C5 of ICarral et all (1 199(H) and 4 - 8 of iBeswick etakl (120031) . 
It is, however, not possible to associate the other two sources, F 
an d G, with any of the sources reported by ICarral et alj (Jl990) 
or lBeswick eTatl d2003l) . 

For all maps, with the exception of the MERLIN L-Band 
one, the parameters (peak coordinates and peak and integrated 
flux densities; Tables [4] and [5]) were derived fitting a 2- 
dimensional Gaussian (AIPS task: JMFIT). In the L-Band map, 
positions were measured using the AIPS task IMSTAT, which 
reads a portion of an image and provides the location and value 
of the maximum brightness in that portion. Integrated flux densi- 
ties were obtained using the task BLSUM, which sums the pix- 
els' flux densities over blotch polygonal regions. The blotch ar- 
eas over which to integrate where determined according to the 
average dimensions of the source (as derived from the other 
maps) in order to avoid confusion with the extended background 
present in the L-Band map. The flux density uncertainties re- 
ported in the two tables (0 & [5]l are provided directly by the 
tasks used and range between a few percent up to 20 percent of 
the values measured. 

At 6 cm, the longest baseline, B max , of MERLIN (215 km) 
provides a resolution of ~ 0.05" that is a factor of 6 better than 
that of the VLA in its most extended configuration. Therefore, 
MERLIN maps can resolve a source seen by the VLA as one in 
two or more components (e.g. our sources A and C). In addition, 
since the shortest baseline, B m \ n , of MERLIN is 9 km, while that 
of the VLA A-array is 0.68 km, the two instruments have quite 
a different response to extended structure (see also Table [2j: at 
6 cm, the MERLIN interferometer is "blind" to emission more 
extended than 1" while the VLA can map emission on a 10"- 
15" scale. This latter problem, which is obvious when compar- 
ing our 5.0GHz MERLIN map with the others (Fig.©, will af- 
fect the computation of the peak flux densities and integrated 
fluxes. The background of the MERLIN L-Band and the VLA 
maps is resolved-out in the C-Band where, rather, a bowl-like 
depression is created , which reduces th e values of the peak flux 
densities (see e.g., IWills et aT1ll998l) . Furthermore, some ex- 
tended emission in the larger of the compact sources may also 
be missed in the calculation of the integrated fluxes. For this rea- 
son, we have decided not to use the 5.0 GHz MERLIN data in 
our computation of the spectral indices. In addition, due to the 
lower resolution of the L-Band MERLIN images some sources 
are blended, and hence, we have also not used this image for the 
spectral index estimates. Nevertheless, it is useful to show the 
MERLIN maps in order to compare the global morphology of 
the radio emission from NGC 520 at different frequencies and 
to potentially spot new sources (e.g. new radio supernovae) that 
might serendipitously appear in galaxies with high star forma- 
tion rates. 

Spectral indices were then computed from the VLA maps at 
8.5 and 14.9 GHz following the convention S oc V . The resolu- 
tion of the two maps is almost identical, and therefore, a convo- 
lution to the same beam was not necessary. The last two columns 
of Table [5] show the computed spectral indices using both peaks 



P. Castangia et al.: New H9O masers in Seyfert and FIR bright galaxies. II. 



7 




28.5 



01 24 35.10 35.05 35.00 34.95 34.90 34.85 34.80 34.75 34.70 

RIGHT ASCENSION (J2000) 



03 47 31.0 
30.5 



O 



O 

LU 

a 



30.0 
29.5 
29.0 h 



I 

C-band 



5.0 GHz 



, o o 



o 



01 24 35.05 35.00 34.95 34.90 34.85 34.80 34.75 34.70 

RIGHT ASCENSION (J2000) 




28.5 



01 24 35.05 35.00 34.95 34.90 34.85 34.80 34.75 34.70 

RIGHT ASCENSION (J2000) 




01 24 35.10 35.05 35.00 34.95 34.90 34.85 34.80 34.75 34.70 

RIGHT ASCENSION (J2000) 

Fig. 4. Radio continuum contour maps of the nuclear region of NGC 520 with labels indicating the compact sources A-G. From top 
to bottom: L-band, MERLIN, contour levels are -3, 3, 6, 9, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48 xO.l mJybeam" 1 ; 
C-band, MERLIN, contour levels are -3, 3, 4, 5, 6 x 0.1 mJy beam" 1 ; X-band, VLA, contour levels are -1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 
10, 11, 12, 13, 14, 15, 16 x 0.1 mJy beam" 1 ; U-band, VLA, contour levels are-3 3456789x 0.1 mJy beam" 1 . The cross marks 
the position of continuum source K detected at 22 GHz (Fig.0. 
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Table 4. Properties of the nuclear continuum sources in NGC 520 observed with MERLIN at 1.4 and 5 GHz. 



Label 


a 


6 


y = 


1.4GHz 


v = 5.0GHz 




(J2000) 


(J2000) 


Peak flux 


Total flux 


Peak flux 


Total flux 




01''24"' 


03°47' 


(mjy/n b ) 


(mJy) 


(mjy/£2 b ) 


(mJy) 


A 


34 s :78 


29 '.'8 


3.5±0.1 


4.7±0.3 


0.56+0.08 


1.0+0.2 




34182 


29 '.'8 






0.55+0.07 


0.55+0.07 


B 


34185 


29 '.'9 


4.5+0.1 


4.5+0.1 


0.38+0.07 


0.38+0.07 




34:86 


29 '.'9 






0.35+0.07 


0.35+0.07 


C 


34:90 


29 '.'8 


4.8+0.1 


5.4±0.4 


0.69+0.08 


2.3+0.3 


C 


34:91 


29 '.'9 






0.54+0.07 


0.54+0.07 


D 


34:92 


29 '.'8 


<0.3 




<0.2 




E 


34193 


29 '.'8 


4.2±0.1 


4.2±0.1 


<0.2 




F 


34 s :97 


29 '.'9 


3.4±0.1 


3.4±0.1 


0.62+0.07 


0.62+0.07 




34:98 


29 '.'8 






0.62+0.07 


0.62+0.07 


G 


35 s :oi 


29 '.'7 


3.3+0.1 


3.3±0.1 


0.44+0.07 


0.44+0.07 



Table 5. Properties of the nuclear continuum sources in NGC 520 observed with the VLA at 8 and 15 GHz. 



Label 


a 

(J2000) 
01''24"' 


6 

(J2000) 
03°47' 


y = i 
Peak flux 
(mJy/fl b ) 


.5GHz 
Total flux 
(mJy) 


v = 14.9GHz 
Peak flux Total flux 
(mJy/n b ) (mJy) 


Or* 
15,p 

Spectral index 
(from peaks) 


<, 
Spectral index 
(from integrated fluxes) 


A 


34:78 


29 '.'8 


1.18+0.03 


5.7+0.2 


0.7+0.1 


2.2+0.1 


-0.8+0.2 


-1.5+0.4 


B 


34185 


29 '.'9 


0.97+0.03 


0.97+0.03 


0.6+0.1 


0.6+0.1 


-0.8+0.3 


-0.8+0.3 




34:86 


29 '.'9 






0.6+0.1 


0.6+0.1 






C 


34190 


29 '.'8 


1.27+0.03 


4.2+0.1 


0.9+0.1 


2.5+0.4 


-0.5+0.2 


-0.8+0.3 


D 


34192 


29 '.'8 


1.29+0.03 


1.29+0.1 


0.8+0.1 


0.8+0.1 


-0.8+0.3 


-0.8+0.3 


E 


34193 


29 '.'8 


<0.1 




0.6+0.1 


0.6+0.1 






F 


34197 


29 '.'9 


0.90+0.03 


5.5+0.2 


<0.3 








G 


35101 


29 '.'7 


<0.1 




<0.3 









and integrated fluxes. The errors on the spectral index a were 
derived using the standard equation for error propagation, 

igv 2 -igv! \l s vi i \ s V2 r 

where v, S v , and AS V are the frequency, flux density, and error 
of the flux density, respectively. The subscripts '1' and '2' refer 
to frequencies 8.5 and 14.9 GHz, correspondingly. 

4. Discussion 

4. 1 . The FIR sample 

The FIR sample of intermediate bright galaxies 
(30<Sioo J um <50Jy) comprises 40 spirals and a galaxy 
pair (NGC 520). Among them, 14 are starburst galaxies and 
only 10 are classified as Seyferts. A s for the FIR sample of 
brighter galaxies (S 100 > 50 Jy; see iHPTl) . there are very few 
galaxies with FIR luminosities greater than 10 n L0 in our sam- 
ple. There are only four Luminous Infrared Galaxies (LIRGs, 
with L FIR >1O U L ), NGC 1614, NGC5135, NGC7469, and 
NGC 7771 and one Ultraluminous Infrared Galaxy (ULIRG, 
Lfir >10 12 L o ), UGC8058. 

Adding the result of our 'intermediate-range' sample, i.e. 
two new maser detections (N GC 61 3 and NGC 520) out of 
41 galaxies searched, to that of IHPTl the cumulative detection 
rate above a 100 fim IRAS Point Source Catalog flux of 30 
and 40 Jy for the parent galaxy becomes 14% and 20%, respec- 
tively (Fig. [5]). When considering the detection rate related to 
the specific FIR flux density interval (30-50 Jy) of our sam- 
ple, we obtain 2/41 or 5%. The two detected galaxies, NGC 613 
and NGC 520, are the second and third FIR-brightest galaxies of 
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Fig. 5. H2O maser cumulative detection rate above a certain 
100 jjm IRAS Point Source Catalog flux density f or the parent 
galaxy. The figure is a modified version of Fig. 13 of lHPTl Filled 
squares represent our values derived using the "intermediate- 
luminosity" sample. 



the sample. The FIR-brightest galaxy of the sample, UGC 2866, 
thus might also be detectable if frequent monitoring is performed 
(nota bene the high degree of variability we see in NGC 520!). 

[HPT! found a rough proportionality between Lfir and Lh,o 
for extragalactic water masers. When plotting the value for 
NGC 613 into this diagram (Fig. [6]), it follows the relationship 
re presenting the Lfir /Lh,q correlation for Galactic masers found 
bv lJaffe etatl dl98ll) . The point representing NGC 520 is, how- 
ever, located well below the line of the correlation. This system 
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Fig. 6. H2O vs. FIR luminosity plot including our two new 
maser detections, NGC 613 and NGC 520, and the separation 
between the kilomaser and megamaser regimes (dashed line). 
Open squares represent most o f the 22 GHz extragalactic H2O 
masers detected so far (fHPTl) . The asterisks indicate the re- 
maining 10 maser sources belonging to the FIR subsample with 
S ioo/im >50Jy (see iHPTf) . The so lid line marks the correlation 
for Galactic star-forming regions (Jaff eet alj|1981l) . The upper 
and lower dot-dashed diagonal lines show the FIR/H2O lumi- 
nosity relations for the mega- and kilomasers , respe ctively, with 
the slope being fixed to that of the Jaff e et ail dl98ll) correlation. 
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Fig. 7. Positions of the main centers of activity in the nuclear 
region of NGC 520. The plus signs mark the radio continuum 
sources C, D, and K detected in our VLA maps, while the five- 
pointed star indicates the position of the H2O maser (this paper). 
The triangl e identifies the kinema tic center derived from CO and 
Hi studies dYun & H ibbard 2001). The circle indic a tes the lumi- 
nous X-ray source (labeled 12) detected by iReadl ((2005) using 
Chandra. The sizes of the marks denote absolute positional ac- 
curacies. 



has either an excess in Lfir w.r.t. Lh 2 o or its maser luminosity is 
particularly low. Most likely, the bulk of the 22 GHz H2O emis- 
sion remains undetected because of limitations in sensitivity (see 
the Sect. l3.2.TT i. When taking a closer look at the Lfir/Lh 2 o plot 
and considering megamaser and kilomaser sources separately 
(above and below the dashed Lh 2 o = 10L© dividing line), the 
two groups seem to follow independent correlations. While H2O 
kilomasers are "subl uminous" for a g iven Lfir w.r.t. Galactic 
star-forming regions dJaffe et alj|1981h . the megamasers appear 
to be slightly "superluminous". The two samples diverge by 
~ 1.5 orders of magnitude in maser luminosity for a given FIR 
luminosity, emphasizing their different origin. 

4.2. NGC 520 

4.2.1. Where is the nucleus? 

Our measurements place the kilomaser in NGC 520 in the inner- 
most nuclear region of the galaxy. In order to associate the maser 
emission with either enhanced star formation or with nuclear ac- 
tivity, we have first to locate and identify the galaxy's true nu- 
cleus and check whether it exhibits (low-luminosity) AGN ac- 
tivity. So far, no clear evidence for the presence of an AGN 
in NGC 520 has been reported. It is impossible to determine 
the spectral type of the main nucleus on the basis of its op- 
tical spectrum becau se of a lack of diagnostic emission lines 
dVeilleux et al.|[l995h . This is likely due to the large extinction 
suffered by the main nucleus (see Sect. I3.2I ). However, an in- 
dication of absent (or weak) nuclear activit y is the presence o f 
PAH features in the mid-infrared spectrum (Bra ndl et alj |2006). 
which should be di luted by the strong infrared continuum of an 
AGN (lKrolikll200Tb . No Fe-K line seems to be present in the X- 
ray spectrum but an ULX source (source 12, see Sect. 13.2b with 
a "power law spectrum" (A^h=2.3 xlO 22 cm -2 ) is detected close 



to the main nucleus (Read 2005). Fig. [7J shows a plot in which 
the positions for the main centers of activity are marked. The X- 
ray source 12 lies c l osest to radio source C (see Sect. l3.2.3l l and, 
according to IReadl d2005l) . is "undoubtedly associated with the 
primary nucleus of NGC 520". However, due to the uncertainty 
of Chandra's absolute astrometry the error circle (« 0'.'7) also 
encompasses sources D and K. Which is, then, the true nucleus 
of the galaxy? 

Source K is close to sources C (further resolved into a dou- 
ble source in the MERLIN C-Band image) and D. Assuming the 
conservative accuracy of O'.'l ( = 15 pc at a distance of 30 Mpc; 
see Fig|4j for the absolute positions in all the VLA maps, source 
K cannot be associated with source C, but an association with 
D is possible (the angular distance is ~0'.'07; see Fig. |7). The 
slightly -resolved source C has a spectral index between 8.5 and 
14.9 GHz, a* 5 , of -0.5 + 0.2 (from the peak flux densities) and 
of -0.8 ± 0.3 (from the integrated fluxes; see Table Q. If the 
slope of the spectrum still holds at higher frequencies, it is pos- 
sible to derive the expected integrated intensities for source C at 
22.5 GHz, 1 .8 +0.4 mJy. In order to make the integrated fluxes for 
source C and K consistent, a spectral index a\ 5 2 of « -4 would be 
necessary, which is impossible to justify with standard physical 
models. For source D, the spectral index is a\ s = -0.8+0.3, pro- 
viding an expected value for the 22.5 GHz integrated intensity of 
0.6+0.1 mJy. This value agrees within the uncertainties with that 
of source K (0.5+0. lmJy; see Sect. 13.2.21 . making the associa- 
tion between the two sources, D and K, the most convincing op- 
tion. Given this interpretation, it remains quite surprising that no 
other source, and especially source C, is detected at 22.5 GHz. 
The non-detection of source C must be a consequence of its an- 
gular extent, confirmed also by the MERLIN C-Band map. The 
other sources are likely too weak to be detected in our highest 
frequency map due to their non-thermal nature. 
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What, then, is the nature of source K (or D)? Its size is 
smaller than O'.'l (= 15 pc) and its spectrum is characterized by 
a steep radio spectral index, which is consistent with optically 
thin synchrotron emission from a RSN or a young SNR, similar 
to the situation found i n other starburst galaxies like M 82 (e. g., 
iMcDonald et"ai]|2002l) and NGC 2146 dTarchi et alj |2000). The 
luminosity of source K (or D) at 5 GHz (extrapolating the flux 
density from the 8.5 GHz value using a spectral index of -0.8; 
Tabled is m 2x 10 20 W Hz -1 , which is greater than the luminos- 
ity of the strongest non-thermal compact source found in M 82, 
but comparable with the luminosities of SNRs found in other 
merger systems, like Arp299 and Arp220 (Fig. [8). However, 
when plotting the surface brightness of source K (~8 x 10~ 17 
W Hz -1 irT 2 sr -1 ) as a function of its diameter (15 pc) in the E-D 
diagram for the SNRs in M 82, NGC 2146, Arp 299 and Arp 220 
(Fig. |5), the value is incompatible with those obtained for the 
other sources (with the exception of the SNRs in Arp 220). 
Being spatially unresolved, the derived source size is an upper 
limit, and hence, the position of source K might move toward 
the upper-left corner of the plot, where the strongest sources in 
M 82, NGC 2146, and Arp 299 lie. Nevertheless, the anomalous 
position of source K in the E-D plot might also suggest that it is 
a low luminosity AGN (LLAGN). Indeed, radio observations of 
Seyfert and LINER gala xies show that LLAGN may have steep 
radio spectral indices (Mor ganti et al.|[l999t iRrips et al.ll2007l) 
and lu minosities comparable to that of source K (Morganti et al. 
1999). In summary, the parameters derived for source K (or 
D), viz. spectral indices, luminosities and diameter, are consis- 
tent both with those of SNRs and LLAGNs. Multi-epoch and/or 
multi-frequency VLBI observations would be necessary to select 
between the two options. 

Source C has also a steep spectral index, a large luminosity 
(2.3 x 10 20 WHz _1 at 5.0GHz) and a position in the S-D dia- 
gram which does not match those of the other sources (Fig. [9]). 
However, one of the key points that determine its position in 
the aforementioned diagram is its size, which we estimate to be 
of order 0'.'3 {=40 pc). This size is too large to justify a possi- 
ble LLAGN nature, and hence, we think that source C is a par- 
ticularly active star-forming region composed of a collection of 
small and young SNRs/RSNe that we do not fully resolve. This 
is also suggested by our C-Band MERLIN map that shows the 
presence of at least one further emission spot (C*, Table|4]i close 
to what we identify as source C. This spot has a size ^0'.'2 (= 
30 pc), a 5 GHz luminosity of * 5.5 x 10 19 WHz" 1 , and a sur- 
face brightness <: 5.7 x 10~ 18 Wm 2 Hz _1 sr -1 . These properties 
are more consistent with those of SNRs (see Figs. [8] and [9), sup- 
porting our interpretation. Also, the kinematic center, although 
close-by, is not coincident with source C favoring a different 
view from that suggested by Read] J2005I) . Given this, it still re- 
mains unclear where the nucleus of the galaxy is located, but if 
a LLAGN exists in NGC 520, it is likely associated with source 
D. 

4.2.2. The H 2 maser in NGC 520: star formation or nuclear 
activity? 

The association of the maser with source K (or D) is based 
on the assumption that the uncertainty in the absolute positions 
is O'.'l in our 22 GHz VLA maps (see Sect. 13.2.21 ). If a less 
conservative value is assumed, line and continuum do not co- 
incide anymore. Therefore, the maser might not be directly as- 
sociated with any of the continuum sources. Nevertheless, the 
small projected distance between the maser spot and source K 
suggests that the H2O emission is in some way related to this 




log(P 5GHz W/Hz) 

Fig. 8. 5 GHz luminosity distribution of the non thermal sources 
C, C, and D (or K) (for this source the flux density has been 
extrapolated from the 8 GHz VLA map using the spectral in- 
dex reported in Table. [5]) in NGC 520 (thi s paper), compared 
with the SNRs in M 82 JmuxIow et al J 1994b. NG C 2146 (Tarchi 
2001, Ph.D. The sis), Arp 299 dNeff et al.1 (20041) . and Arp 220 
dParra et alj|2007h . 
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Fig. 9 . 2-D relation for SNRs in M 82 (crosses; iMuxlow et al.l 
1994), NGC 2146 (squares; Tarchi 2001, Ph.D. Thesis and 
Surci s, 2004, Ph.D. Th esis), Arp 299 (t riangles; iNeffetalJ 
l2004b . Arp 220 (asterisks: IParra et al]l20"07b . and the sources C, 
C*, and D (or K) (diamonds) in NGC 520 (this paper). The arrow 
shows the direction into which source D would be moved, if its 
size would be <15 pc. 



source. The characteristics of the maser line, i.e. narrowness and 
strong time variability (70% decrease of the flux density within 
24 days) are more consistent with star formation than with nu- 
clear activity. Narrow flaring components are common in wa- 
ter maser emission from Galactic star-forming regions, although 
only the outbursts in W49N reached a luminosi t y comparable 
to that of NGC 520 dLiriestrom & Gwinnl l2000t iHonma et al.l 
12004 . Narrow extragalactic H2O maser lines ha ve been de- 
tecte d, so far, also in IC 10 dHenkel et al.1 1 1986b and IC342 
dTarchi et al.ll20 02b). Both these kilomasers are definitely asso- 
ciated with star-forming regions and have shown a great deal 
of time variability. The flare of the -324km s _1 component in 
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IC 10 reached ~4L Q or 120Jy dBaan &Haschicklll994i) onl y 
to decay to sub-Jy levels more recently (Brunth aler et al.l l2007). 
Even more similar to the H2O maser in N GC 520 is the maser 
in IC 342 reported bv lTarchi etall d2002bl) . Although less lumi- 
nous (Lh 2 o = 10~ 2 Lq), the feature was n arrow (~ 0.6 km s' 1 ) 
and varied by a factor of 3 within 5 days (iTarchi et al.ll20 02b). 
Moreover, for the maser in IC 342, a velocity shift by 1 km s _1 
was observed within two weeks. The shift and the rapid decay in 
flux was explained in terms of a chance alignment of two dense 
molecular clouds. 

A similar scenario is also compatible with the maser in 
NGC 520. When comparing our highest velocity resolution spec- 
tra taken with Effelsberg on January 2005 with that taken with 
the VLA a month later, a small velocity shift seems present. A 
Gaussian fit to the maser lines indicates velocities of 2275.42 
± 0.08 and 2275.15 + 0.03 for the Effelsberg and the VLA, re- 
spectively (Fig. [2]». The velocity shift is 0.27+0.09 km s _1 , i.e. 
it may be significant. We may assume that both the flare and the 
shift (if re al) have a kinematic origin, following the scenario pro- 
posed by Tarc hi et ail ((2002b). Assuming a negligible distance 
between the two clouds w.r.t. the distance of the galaxy and a 
relative velocity in the plane of the sky of ~ 100 km s _1 , this im- 
plies a cloud velocity gradient of up to 0.23 km s~'/AU during 
our monitoring time (24 days). Alternatively, the maser feature 
may be explained by a change in the relativ e intensities of the 
three main hyperfine components (see e.g., Fiebig &" Guested 
Il989l) or by a variation in the relative intensity of different emit- 
ting regio ns leading t o a swapping-over of velocity peaks (see 
e.g., IWu et al.lll999l) . For maser flux density variations (dis- 
tinct from flares) due to refractive interstellar scintillations, time 
scales range from about one month (slightly long er than the vari- 
ability time scale in NGC 520) to few months ( Simonett i et alj 
1 1 993b . Diffractive scintillation may instead allow for extremely 
rapid variations (as short as a few minutes) but ha s been reported, 
so far , in only one case, viz. the Circinus Galaxy (Greenh ill et alj 
119971) . 

Although the attributes of the maser (low luminosity, nar- 
rowness, and strong variability) hints at an association with star 
formation, we cannot rule out the possibility that we have found 
a weak nuclear kilomaser, as the ones in M 5 1 and (possibly) 
NGC 4051, perhaps related to a nuclear jet that is too weak to 
be seen in our radio maps. However, no high-velocity lines blue- 
or redshifted with respect to the systemic velocity are detected 
in NGC 520 and this does not support the possibility that the 
maser is tracing an accretion disk. In addition, the linewidth of 
the maser feature in NGC 520 (~ 0.6 km s ) is very small com- 
pared to the typical linewid th of H2O maser s assoc iated with 
nucle ar jets (e.g., Mrk 348. IFeck et alj l200l M51. lHagiwaral 
2007). Most of all, no convincing evidence of an AGN of what- 
ever intensity comes from any study performed so far, making 
the latter alternative unlikely. However, also a few megamaser 
galaxies, seen edge on and thus being extremely obscured (like 
NGC 520), do not show th e presence of an AGN JGreenhill et alj 
l2003btlBraatz et alj|2004 . 

5. Conclusions 

We have observed a sample of 41 galaxies with 100 /jm IRAS 
point source flux density between 30 and 50 Jy and declination 
> -30°, to search for 22 GHz H2O maser emission. The survey 
yielded two detections, a megamaser with isotropic luminosity, 
L H ,o, of « 35 Lp in th e Seyfert/HII galaxy NGC 613 (see also 
iKondratko et alj|2006l) and a kilomaser with Lh,o~ 1 L0 in the 
merger system NGC 520. Our results confirm the correlation be- 



tween maser detection rate and FIR flux density discovered by 
iHPTi for the sample of galaxies with S 100/im > 50 Jy. When we 
compare the isotropic luminosity of extragalactic masers with 
the FIR luminosity of their parent galaxy, we find that kilomasers 
are "subluminous" while megamasers are "superluminous" for a 
given Lfir w.r.t Galacti c H2O masers i n mas sive star-forming 
regions. We agree with IKondratko et al.l d2006l) that the maser in 
NGC 613 is possibly associated with the nuclear jet, although the 
fact that the H2O emission is centered on the systemic velocity 
is unusual for a jet-maser and may hint at a different origin. 

The kilomaser in NGC 520 has been observed at subarcsec- 
ond resolution with the VLA. We have found that the H2O emis- 
sion arises from the inner hundred parsecs surrounding the dom- 
inant nuclear region of the merger system. Assuming an abso- 
lute position accuracy of 0'.' 1 (= 15 pc), the maser coincides with 
a continuum source (source K, or D). A spectral index analy- 
sis, as well as its general properties, i.e., its size and luminos- 
ity, indicate a young SNR or a LLAGN. The small linewidth 
(~ 0.6 km s _1 ) and the time variability (70% decrease of the flux 
density within 24 days) of the maser feature suggest a size limit 
of $0.02 pc and a brightness temperature greater than 10 10 K, if 
the fluctuation is not caused by a variable background contin- 
uum source. The H2O maser emission is most likely associated 
with star formation, like the majority of kilomasers. Narrowness 
and variability of the maser feature resemble those of two other 
kilomasers observed in IC 10 and IC 342 that are associated with 
star formation. We consider an association with an AGN a less 
viable option based on the lack of clear traces of AGN activity. 
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